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1 Introduction 



In this paper we consider the microscopic and the macroscopic properties of a class of 
generalized statistics, referred to as A-superstatistics. In our approach, the starting 
point is a certain symmetry principle, characterized by an algebra of creation and 
annihilation operators (CAO's) with Fock type representations. 

The idea behind these investigations is based on a few observations. The first 
one is that any n pairs of Bose CAO's Bf, ■ ■ ■ , Bf generate a representation, 
the Bose representation p#, of the orthosymplectic Lie superalgebra osp(l\2n) = 
B(0, n). The representation independent generators Bf, Bf, ■ ■ ■ , Bf of osp(l\2n), 
which in the Bose representation coincide with the Bose operators ps{B~f) = Bf, 
are para-Bose CAO's pQ and satisfy the relations 

[{BlB]},Bl\ = (e-05 ik B] + (e-r))5 jk Bf, £, 77, e = ± or ± 1. (1.1) 

These triple relations give one possible definition of the Lie superalgebra osp(l\2n) [2], 
and are considered as the defining relations of para-Bose statistics. So ordinary Bose 
statistics corresponds to one particular realization of para-Bose statistics. 

The situation with Fermi statistics and its generalization, para-Fermi statis- 
tics pQ, is similar. So the second observation is that any n pairs of Fermi CAO's 
Fi, F«f , • • • , F„ give a representation, the Fermi representation p F , of the orthogonal 
Lie algebra so(2n + 1) = B n and any n pairs of para- Fermi CAOs Pf, with 

[[^,F^F^ = l( V -e)% k ^-^-e)% k F^ £, V ,e = ±or ± 1, (1.2) 

generate the algebra B n j^j. 

Both algebras B n (every Lie algebra is also a Lie superalgebra) and B(0\n) belong 
to the class B of basic classical Lie superalgebras in the classification of Kac j3]. 
Hence ordinary Bose and Fermi statistics and their generalizations para-Bose and 
para- Fermi statistics could be referred to as B- (super) statistics. 

On the ground of these facts statistics related to the other classes of basic Lie 
superalgebras were introduced : A-, B-, C- and _D-(super)statistics 0. So far only 
A-statistics (corresponding to the Lie algebra sl(n + 1)) was studied in detail from 
the microscopic [5], jSl and macroscopic point of view [7j. In the present paper we 
investigate the properties of A-superstatistics jH], namely the statistics arising from 
the Lie superalgebra sl(l\n) = A(0, n — 1) or gl(l\n). 

The microscopic properties of this statistics related to gl(l\n) lead to the in- 
troduction of Fermi-like operators , and these are relevant for the description of 
certain physical models (see Section 3). In particular, lattice models of strongly 
correlated electron systems use these operators ff (in a certain representation). 
Also, the so-called ideal odd-particle creation and annihilation operators of Klein 
and Marshalek in nuclear shell model theory correspond to a particular realization 
of the operators ff. 

Section 2 is devoted to the microscopic properties of this new class of generalized 
statistics. First we recall the definitions of the Lie superalgebras sl(l\n) and gl(l\n) 
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and their Fock representations. Like for para-Bose and para-Fermi statistics, the 
generators f^, i = 1, . . . , n of sl(l\n) satisfy certain triple relations. These relations 
completely define the Lie superalgebra sl(l\n), just as the triple relations for the 
generators af, i — 1, . . . ,n of A-statistics do for the Lie algebra sl(n + 1) 0. For 
sl(n + 1), this property of the operators af was first observed by Jacobson and 
therefore the generators af are referred as Jacobson generators of sl(n + 1) jE]. By 
analogy we call the generators ff 1 of sl(l\n) also Jacobson generators of sl(l\n). 
The Fock representations of A-superstatistics are constructed in the same way as 
in parastatistics. They are generated by the operators ff~ and labeled by a posi- 
tive integer p = 1,2, ... . Within the corresponding module W(p,n) the generator 
ff~ (respectively ff) can be interpreted as a creation (respectively annihilation) 
operator of a "particle" on the zth orbital (i = 1,2, ... ,n). The Pauli principle 
of A-superstatistics is formulated and in Section 3 it is indicated that in the limit 

p — > oo the representation dependent operators F(p) i = coincide with ordinary 
Fermi CAO's. We complete this section by showing that the Jacobson generators of 
sl(l\n) are implicitly present in certain physical models. 

The following sections are devoted to the macroscopic properties of A-super- 
statistics. In Section 4 we construct explicitly the sl(l\n) grand partition function 
Z(p,n) (GPF), the average number of particles in the system N(p,n) and the av- 
erage number of particles on each orbital 9i under the assumption that the energy 
of each particle on orbital % is €{. All these thermal properties of the system are 
described by means of the elementary symmetric functions. The fact that symmet- 
ric functions appear naturally in the description of (grand) partition functions in 
statistical mechanics got attention recently ^U] • For quantum systems with Bose or 
Fermi statistics, see ^Uj- For quantum systems with A-statistics, the relevant sym- 
metric functions are the so-called complete symmetric functions [Zj. Here in the case 
of A-superstatistics, the relevant symmetric functions are the so-called elementary 
symmetric functions. 

Then we consider two specializations of the general case of Section 4. In Section 
5 the energy and the chemical potential of each orbital are assumed to be the same. 
This is the so-called degenerate case. The thermo dynamical functions simplify, 
and many of these can be expressed as hypergeometric series. Furthermore, these 
functions can be seen as deformations of the corresponding ones in the case of Fermi 
statistics. In Section 6 we investigate a model with equidistant energy levels. Also 
under this specialization, the thermodynamical functions assume a simple form, 
usually in terms of g-generalized or basic hypergeometric functions. 

It should be emphasized that our notion "A-superstatistics" refers to a gener- 
alized quantum statistics based upon the Lie superalgebra of type A, and that it 
should not be confused with the recently introduce phenomenon of "superstatistics" 
as in jTTj. 
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2 Microscopic properties of A-superstatistics 



A-superstatistics is defined in the context of the Lie superalgebra sl(l\n) or gl(l\n). 
A convenient basis of gl(l\n) is given with the Weyl generators e^, with i,j G 
{0, 1, . . . , n}. The grading of gl(l \n) is as follows : the even elements are given by 
eoo and with i, j e {1, . . . , n}; the odd elements are e^ and e^ {i = 1, . . . , n). 
The Lie superalgebra bracket is determined by 

[ey, e w ] = eye w - (-lf^^^e^ = 5 jk e u - (~l) dc ^ dE ^5 u e kj , (2.1) 

where deg(ejj) is (resp. 1) if is even (resp. odd). One can define sl(l\n) as the 
(super)commutator algebra of gl(l\n); its basis consists of all elements ey (z 7^ j) 
and the Cartan elements e 0Q + en (i = 1, . . . , n). 

The Jacobson creation and annihilation operators ff of sl(l\n) are given by 

ft = e *o, /r = e o*> (i = 1, . . . , n). (2.2) 
It is known [Hj that the linear envelope of 

{# {/?, fk}\h j,k = l,...,n;t,r,,e=±} (2.3) 

is indeed the Lie superalgebra sl(l\n). 

A-superstatistics is determined by the relations that hold for the creation and 
annihilation operators. These are : 

{ftJt} = {fr>fi} = °> ( 2 - 4 ) 

[{ft,fi}Jt}= Sjkft-Saft, (2-5) 
i{ftJf}Jk] = -tikfi + tiifk- (2-6) 

It is worth observing that these operators satisfy the compatibility conditions re- 
quired in the context of Wigner quantum systems ^2] ■ 

The Fock representations of A-superstatistics have been classified by Palev jH]. 
These representations are labeled by a positive integer p, the order of statistics. 
Let us denote the Fock representation of order p for sl(l\n) by W(p, n). The space 
W(p, n) is characterized by a vacuum vector |0), such that 

/T|0> = 0, (i = l,...,n), (2.7) 
Uft\0)=p5 lJ \0), [i. j 1 n). (2.8) 

These Fock spaces are finite-dimensional unitary irreducible s/(l|n)-modules. A set 
of basis vectors for the space W(p,n) consists of all vectors 

(ft) 6l (m d2 ---(ft) 9n \0), * 4 e{0,l} (2.9) 

where 

n 

\e\ = J2^<p- (2.10) 

i=l 
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The linear span of all vectors (|2.9jl . without the restriction ([2.10)1 . also forms a 
sl(l\n) module W(p,n). However, if p < n, W(p,n) is not irreducible : it contains a 
maximal invariant sub module, and W(p, n) is the quotient module of W(p, n) with 
respect to this submodule. If p > n, we have that W(p,n) = W(p,n); in this case 
it is clear that the restriction ()2.10j) is superfluous. 

One can define a Hermitian form ( , ) on W(p,n) with the usual Fock space 
technique, by requiring 

(0|0) = 1, (2.11) 
(ftv\w) = (v\f?w), Vv,weW(p,n). (2.12) 

With respect to this form, the different vectors in (|2.9|) are orthogonal, and the 
following vectors form an orthonormal basis of W(p, n) : 



\ P - e) = \p- e lt . . . , e n ) = J ^P^(tf) ei (m 62 ■ ■ ■ (/n + ) e lo), 

^£{0,1}, \9\<p. (2.13) 

Furthermore, the Hermitian conjugate of / 4 is in this module, which is an 
important physical requirement. 

The transformation of the basis ()2.13|) under the action of the creation and 
annihilation operators reads : 

fr \ p; e) = ^(-1)^+-+^ y/ p - \e\ + i \ v - e 1 ,...,e i -i,...,e n ), (2.14) 
f t +\m e) = {i- \p;0 1 ,...,e i + i,..., e n ). (2.15) 

Our Hamiltonian will be an element from (the Cartan subalgebra of) gl(l\n). There- 
fore we first extend W(p, n) to a gl{l\n) module. For this purpose, we set iVj = en 
(i = 0, 1, . . . , n), and have : 

N Q \p-e) = (p-\6\)\p-e), (2.16) 
N l \p-e) = e l \p-6). (2.17) 

We shall be studying macroscopic properties of A-superstatistics for a Hamilto- 
nian of the following form : 

n 

H = J2^. (2.18) 

i=l 

Via creation and annihilation operators, this can be rewritten as 

n / n \ 

H = E e * /r> + — x (p - J2if^ f k }) ■ ( 2 - 19 ) 

1=1 V fc=i y 

Clearly, H\0) = (so the vacuum has zero energy), and 

[HJt]=±eJt- (2-20) 
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Therefore, each (resp. / t ~) can be interpreted as an operator creating (resp. 
annihilating) a particle (or quasiparticle, or excitation) on orbital i (with energy e^). 
Since 

H\ P ]6)= M>> ( 2 - 21 ) 

\p\ 9) is interpreted as a state with &i particles on orbital i. 

The CAO's of sl{l\n) together with the Hamiltonian (j2.18|) generate gl(l\n). 
Then, in view of ()2.19|) . f\,---,fni considered as operators in any W(p, n), generate 
gl(l\n). For this reason we call the CAO's of sl(l\n) also Jacobson generators of 
gl(l\n). 

The Pauli principle for A-superstatistics follows basically from equation (|2.1()|) 
and the form of the vectors (j2.13|) . These relations imply that the system can 
accommodate up to min(p, n), but no more than min(p, n) particles, in such a way 
that every orbital contains at most one particle. 

Note that in the atypical cases (p < n) the Pauli principle introduces an inter- 
action, a statistical interaction, between the orbitals. For instance if the system 
is in the state with 9i = . . . = 9 P = 1 and 9 p+ i = . . . = 9 n = 0, then it cannot 
accommodate more particles, not even on the empty orbitals with i > p. Therefore 
the orbitals are not filled independently. The filling of a given orbital depends on 
fillings of all other orbitals, which is the main feature of exclusion statistics [T5] . 
Observe also that A-superstatistics is closely related to ordinary Fermi statistics : 
the only extra condition comes from the order of statistics p. 



3 Quasi-Fermi creation and annihilation opera- 
tors 

In the present section we discuss some differences and similarities between A-super- 
statistics and ordinary Fermi statistics. First of all, observe that there is an impor- 
tant difference between the Fock spaces W(p,n) with p < n and those with p > n. 
For p < n, the irreducible sl(l\n) representations are atypical, and their dimension 
is given by : 

dimW(p,n) = ^f n \ (3.1) 

fc=0 ^ ' 

For p > n, the representations W(p,n) are all typical, with 

dimW(p,n) = 2 n . (3.2) 

For pi ^ p 2 and Pi,P2 > n , the representations W(pi,n) and W(p2,n) are certainly 
not isomorphic (even though they have the same dimension), since they have a 
different highest weight. 
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Let n be fixed, p > n, and let us define representation-dependent operators 
F(p)f in W(p,n) by 

F(p)f = ^, z = l,2,...,n. (3.3) 

The action of these new creation and annihilation operators on the vectors (J2.13)) 
reads : 



F(p);\p;6) = e i (-i)^^ili + L-M\p-0 1 ,...,O i -l,...,O n ), (3.4) 



F(p)l\p; 9) = (1 - e^-lf^-^Jl - B \p; 9 1; . . . A + 1, . . . , 9 n ). (3.5) 



On the other hand, one can consider a set of n ordinary Fermi creation and annihi- 
lation operators Ff 1 (i — 1, . . . , n), satisfying : 

{Ft,F;} = {F-,Fr} = 0, {Ff, = (3 . 6) 

and its 2 n -dimensional Fock space W(n) with orthonormal basis vectors 

\9) = \9 U . . . , 9 n ) = (F+f ■ ■ ■ (F+) e " |0>, 9, e {0, 1}, (3.7) 

and action 

Fr\9) = 9 i (-l)^+- +e ^\9 1 ,...,9 i -l,...,9 n ), (3.8) 

F+\ ) = (l-9 i )(-l) e ^-+ e ^\9 1 ,...,9 i + l,...,9 n ). (3.9) 

We can now identify the basis vectors \p; 9) of W(p, n) with the vectors \9) of W(n). 
Let us therefore define, for a given positive integer p > n, in W(n) a realization of 
the operators F(p)f, denoted by p(F(p)f), through the action : 



p(F(p)r)\9) = 9 i (-l) e i + - +e >-^ll + —^\9 1 ,...,9 l -l,...,9 n ), (3.10) 



p{F{p)t)\9) = (1 - 9i)(—l) ei+ "' +ei ~ 1 \l 1 - M \9t, . . . , 9, + 1, . . . , 9 n ). (3.11) 

Both p(F(p)f) andi^f are now operators acting in the finite-dimensional Fock space 
W(n), and from their action it follows immediately that 

lim p(F(p)f) = F^ (3.12) 

p— +oo 

For this reason, the operators F(p)f are said to be quasi-Fermi creation and an- 
nihilation operators. For large p-values, they tend to ordinary Fermi creation and 
annihilation operators. 
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It is interesting to consider the anti-commutators of the quasi-Fermi creation and 
annihilation operators, acting in the Fermi Fock space Win). From ([3.10|) - ([3.1ip one 
obtains 



{p(F(p)t),p(F(p)f)}\9) = 
{p(F(p)r),p(F(p)j)}\9) = 

{p(F(p)r), p (F(p)i)}\6} 

{p(F(p)r), p (F(p)+)}\e} 

x 

{p(F(p)-),p(F(p)l)}\6) 



1 + ^)|»>. 

P 



. , 8i - 1, . . . , 9j + 1, . . .), 

\-l) e i+-+ d i9 i {\ - 8,) 



i < J, 



(3.13) 
(3.14) 

(3.15) 



(3.16) 



V 



,9 j + l,...,9 i -l,...), i>j. (3.17) 



Compare again with ()3.6|) to see that the above anti-commutators tend to ordinary 
Fermi anticommutators when p tends to infinity. Also in this sense, the quasi-Fermi 
creation and annihilation operators p{F{p)f) can be considered as "deformations" 
of ordinary Fermi creation and annihilation operators, with the integer p (the order 
of statistics) SIS db deformation parameter. 

A comparison of ()2.14j) - (j2.15j) for p > n with ()3.8)) - (j3.9j) shows that the Jacobson 
generators of sl(l\n) can be expressed as functions of ordinary Fermi creation and 
annihilation operators F^, F^, . . . , F^. Let W(p, n) be a typical representation (so 
with p > n), and identify again its basis vectors \p; 9\, . . . , 8 n ) with the basis vectors 
\9i, . . . , 9 n ) = \9) of the Fermi Fock space W{n). Since F^F~ is a number operator 
for fermions in a state i, 



F l + Fr\8) = 8 l \8} 1 



l,...,n, 



(3.18) 



we can write for (J2.15|) 

^l^) = e J o|^) = (l-^)(-l) ei+ 
The latter can be represented as 
f?\0) 



V 



+ l-Y J F^F- \8 u ...A + l, 



k=i 



■ , 9 n ). 
(3.19) 



^-J2 F k F kFt\0) = F^ 



k=l 



\ 



P 



J2 F kK\°)- (3-20) 



k=l 



Equation ()3.20|) holds for any |#). Therefore 



f'i 



eio 



i 



p-J2 f ? f *- 



n. 



k=i 



(3.21) 
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In a similar way one derives from (|2.14|) 



P-J2 F k F kK~, i = l,...,n. (3.22) 

fe=i 

Evidently fsee l2~TH|) . 

n 

e 00 =p-J2FkFk, (3-23) 

fe=i 

and simple calculations lead to 

/•}'/•;• i,j = l,...,n. (3.24) 

In such a way we have expressed all Weyl generators {e^|z,j = 0, 1, . . . , n} of gl(l\n) 
via n pairs of Fermi operators. Let us mention that the "ferminization" ()3.21|) - (j3.22|) 
of the Jacobson generators of sl(l\n) is not new. It is known as the Holstein- 
Primakoff realization of sl(l\n) [T2|. The use of boson or fermion operators to 
construct representations of Lie algebras or Lie superalgebras has a long history. 
For Lie superalgebras, see fHJ IS] for a realization of the algebra in terms of Bose 
and Fermi operators and a construction of the corresponding representations. 

Finally, we wish to mention that the Jacobson generators of gl(l\ri) and the 
considered Fock representations W(p, n) (or equivalently, the quasi- Fermi operators) 
are implicitly present in certain physical models. 

Examples from condensed matter physics include mainly models related in one or 
another way to high-temperature superconductivity. We have in mind those lattice 
models of strongly correlated electron systems where (the electronic part of) the 
Hamiltonian is expressed in terms of Hubbard operators (X-operators) as for 
instance in JH] or in • hi such models, each Hubbard operator is labeled by three 
indices, X^, where A refers to the lattice site and i, j — 0, 1, . . . , N, if the (combined 
spin-flavor) degrees of freedom (the number of the orbitals) of the electrons at each 
fixed site are N. For any site the operators X A and X A l , i — 1, . . . , N, are said to 
be fermion-like generators (or of odd degree), whereas X A °, X A , i,j — l,...,N are 
boson-like generators (or of even degree). The X-operators obey the relations 

[X'l X k B % = 5 AB {5 ik Xl ± 6 u X* j ), i,j,k,l — 0,..., N, (3.25) 

with the upper signs if both X-operators in the left hand side are fermion-like and 
with lower signs in all other cases. Clearly (J3.25|) can also be written as 

[X% X*f\ = 5 AB {5 ]k Xl - (-l) d ^ x ^ d ^ xk ^5 u X k A >), (3.26) 

which indicates, that for a fixed value of the lattice site A the Hubbard operators 
are the Weyl generators of gl(l\N) = gl(l\N) A [20] > cfr. (J2.1j) . Since moreover for 
A 7^ B the X-operators supercommute, [X^,X|j] = 0, the conclusion is that all 



fi 



eoi 



\ 
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Hubbard operators constitute a basis in the algebra C which is a direct sum of all 
gl(l\N) A , i.e., 

C = Qgl(l\N) A . (3.27) 

A 

Each local state space per site A (we suppress the site index A whenever possible) 
has a basis consisting of all vectors 

|n ;ni, ...,n N ), n 1 , . . . ,n N e {0,1}, (3.28) 

subject to the additional constraint n = p — J2k=i n ki with p fixed (in [TH], p = N/2; 
and in [THj, p = 1). Since n is required to be a non-negative integer, only such sets 
of numbers n\, . . . ,n^, are admitted for which 

m + . . . + n N < p. (3.29) 

The physical interpretation of the state (|3.28j) is that it corresponds to a config- 
uration with ri\ electrons on the first orbital, n 2 electrons on the second orbital, 
etc. The action of the X-operators on the states ()3.28|) reads (throughout below 

i, 3 ,k = \,...,N) m ■ 

X k A k \n ;ni, ...,n k ,.. . ,n N ) = n k \n ; n 1 ,...,n k ,.. .,n N ), (3.30) 
X 3 A |n ;ni, . . . ,rij, . . . ,n k , . . . ,n N ) = 

(-l)%+-+"*-i|n ; n 1 ,...,n j + l,...,n k -l,...,n N ), j ^ k, (3.31) 
X A \n Q ;ni,...,n k ,...,n N ) = 

v ^TT(_i)«i+-+^-i | no + i ; nu . . . , n k - 1, . . . , n N ), (3.32) 

^A°K;Wi, ...,n k ,.. .,n N ) = 

^ / _ ( _ 1)ni+ ... +nfc _ i ^ _ 1; ni? _ _ _ ) nfc + 1; _ _ _ j njv ^ (3.33) 



In eqs. (I3.30|) - (j3.33|) . the convention is that the vectors on the right hand sides with 
unacceptable arguments should be identified with zero. 

Equations ()3.32|) and ()3.33|) clearly indicate that the operator (resp. X A k ) 
creates (resp. annihilates) an electron on the A;-th orbital of site A. However, 
these operators are not Fermi creation and annihilation operators in the strict sense 
because {X A \X A } ^ 

Comparison with the formulas from Section 2 now leads to the identification of 
the Hubbard model Hilbert space with the representation W(p,n), where 

N = n, (n l ,...,n N ) = (Q l ,...,6 n ), n = p - \9\. (3.34) 

Furthermore, the Hubbard operators (on a fixed site A) are expressed in terms of 
the Jacobson creation and annihilation operators f A of gl(l\n) : 

Xf = 4", Xf = /+ X A k = {//, /-}, (j < k). (3.35) 

So the conclusion is that the operators X A ° and X A k , creating and annihilating 
electrons at site A, are not Fermi operators. They are creation and annihilation 
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operators of the Lie superalgebra gl(l\n). The statistics of the electrons or of any 
other (quasi)particles described with these operators is not Fermi statistics, it is 
A-superstatistics of order p. 

In order to quote an example from nuclear physics, note that the p — 1 Jacobson 
generators of sl(l\n) together with N satisfy the relations : 

0, (3.36) 

(3.37) 

0, (3.38) 

(3.39) 

In nuclear shell model theory the operators with the above properties are called 
ideal odd-particle (IOP) creation and annihilation operators |21j . Okubo [22| refers 
to the algebra of IOP operators as to the Marshalek algebra. The IOP operators 
play a relevant role for the description of properties of odd nuclei in the frame of 
the nuclear shell model (see the review article [23] and references therein). 



f 7 f 7 = f+ ft = 

ST ft = 

N ft = ffN = 
iV 2 = N . 



4 Macroscopic properties of A-superstatistics 

To describe the macroscopic properties of A-superstatistics for the Hamiltonian H 
with orbitals % (i — 1, . . . , n), see (|2.18|) . it is for us irrelevant whether the different 
orbitals correspond to different particles, to different energy levels of particles, or to 
different internal states of the particles. The only assumption is that they satisfy the 
Pauli principle for A-superstatistics, which follows from the Fock space construction. 

As usually, we assume that the system is in a thermal and diffusive contact and 
in a thermal and diffusive equilibrium with a much bigger reservoir. We denote by 
t its (fundamental) temperature, by the chemical potential and by the energy 
for the particles on orbital i. 

The probability V(p, n; 9) for the system to be in a (quantum) state 9 = (9\, . . . , 9 n ) 
with \9\ = 9\ + ■ • • + 9 n particles and energy E = 9\t\ + • • • + 9 n e n is given by the 
expression 

^) = ^^{t^)- ("J 

The numerator in this expression is the Gibbs factor of the system in the state 9, 
and Z(p,n) is the grand partition function (GPF) of the system. In the case of 
Bose or Fermi statistics or their generalizations (Green parastatistics |T]), the GPF 
is simply a product of the GPFs of all orbitals. This is due to the fact that for 
those statistics the different orbitals can be considered as independent subsystems : 
the filling of each orbital is completely independent of the number of particles that 
have already been accommodated on the other orbitals. Here due to the new Pauli 
principle this is no longer the case if p < n. Therefore we have to compute directly 
the GPF for the whole system. The latter as usual is the sum of the Gibbs factors 
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over all possible states of the system. So we have : 

Z(p,n)= fexp(^^i)) 9i -.-(exp(^^)) . (4.2) 



o<e 1 +---+e„< P 

^£{0,1} 



In terms of the notation 



Xi = exp(— — — ), i = l,...,n, (4.3) 
r 



we have 



min(p,n) 

Z(p,n) = Yl 4 1 4 2 ---x 9 n= J2 T, Xi4 2 '--x 9 n- (4-4) 

0<9]H ye n <v k=0 9i + --- + Sn=fe 

e,G{0,i} ^£{0,1} 

It follows that Z(p,n) = Z(n,n) if p > n; therefore, we shall from now on assume 
that p < n, thus covering all possible cases. Since all macroscopic properties are 
encapsulated in the grand partition function, this observation also implies that A- 
superstatistics with p > n has the same macroscopic properties as ordinary Fermi 
statistics. So the case p = n — 1 can be considered as the smallest deviation from 
Fermi statistics, as far as the macroscopic properties are concerned. In the following, 
we shall sometimes pay special attention to the p = n — 1 case, and compare it with 
the properties of Fermi statistics. 

In the present context, it is appropriate to introduce the elementary symmetric 
functions ek(x\, . . . , x n ), k = 0, 1, . . .. The k-th elementary symmetric function [21] 
is the sum over all products of k distinct variables Xj, so that e (xi, . . . , x n ) = 1 and 



6 k ("El ; • • • ; ^ -^i2 ' ' ' ^ik (4' 5) 

ii<«2<---<«fc 

xl'xl 2 ■■■x 9 n n . (4.6) 



e 1 + ---+e n =k 

«i€{0,l} 



For instance, e 1 (x 1 ,x 2 ,x 3 ) = x 1 + x 2 + x 3 , e 2 (x 1 ,x 2l x 3 ) = x x x 2 + x x x 3 + x 2 x 3 , 
e 3 (xi,x 2 , x 3 ) = Xix 2 x 3 , and e^{xi, x 2 , x 3 ) = for k > 3. The generating function 
for the e/c is given by j21] 

n 

ek(zi, • • • > %n)t k = (1 + xit) ■■■{! + x n t). (4.7) 

In terms of the elementary symmetric functions, one finds 

p 

Z(p, n) = >^ e k (xi, ...,x n ). (4.8) 

k=0 
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This sum does not simplify if p < n; for p = n, f)4.7|) yields 

Z(n, n) = (1 + Xl )(l + x 2 ) ■ ■ ■ (1 + x n ). (4.9) 

Also from ()4.7|) . it follows that one can give the following description : for p < n, 
Z(p, n) consists of those terms of Z(n, n) that have total degree less than or equal 
to p. 

Let us now consider some other thermodynamic quantities. The probability 
V(p, n; 9) for the system to be in the state 6 = (9 1 , . . . , 6 n ) with 16*1 particles reads 

V(p,n;6) = X ^ X ; (4.10) 
Z(p, n) 

Therefore, the average number of particles in the system is 

N(p,n)= £ \0\V( P ,n;9)= £ \9\ 1 2 " . (4.11) 

0<6 1 +—+e„<p 0<9 1 + ---+B n <p ' 

»i6{o,i} Sie{o,i} 

This can be rewritten as 

n 

N(p,n) = J2x t d x An(Z(p,n)). (4.12) 

i=l 

In terms of symmetric functions, the numerator of (|4.11|) reads 

v P 

E\Q\^,0l 02 n — \ U \ „6l 02 r,fin — \ hp, I rr v \ 

— / j ™ / x l x 2 x n ~ / j fcc fcl A l; • • • ; 

0<9 1 +-+9n<P /c=0 ei + ---+e n =fc fc=0 

^€{0,1} ^£{0,1} 

(4.13) 

so we can write 

N(p, n) = ^M*!, ■■■,*») (4 w) 

or equivalently, 

jV>, n) = p - ^(P-fch^i.-^) , (4 . 15) 

From this last equation it is clear that the average number of particles in the system 
is indeed less than p. Formula f)4. 14j) simplifies when p = n. Indeed, from ()4.7|) we 
have 

n d f " 

ke k (x!, x n )t k = t—[^2 e k{xi, %n)t h 

fc=0 \k=0 

q / n \ n n 

= ^ ri( i + ^) = 1 e rrirt n^ 1 + ^ ( 4 - 16 ) 

\i=l / r=l r i=l 
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and so 

n 

r=l r 

Next, we shall determine the equilibrium distribution of the particles on the 
orbitals. First of all, consider the orbital n. Either there are no particles on this 
orbital, or else there is just one particle. Denote by V(p, n; 9 n = 0), resp. V(p, n; 9 n = 
1), the probability that there are no particles present (resp. that there is one particle 
present) on orbital n. From the sum of the corresponding Gibbs factors, one finds : 

P(P , „;»„ = (,) = Z 'P-»)I T ° = XXo^. ■■■■*-.) (4 . 18) 

Z(p,n) ^ =0 e fc (xi, . . . ,x„_i,x n ) 

^-1 



V(p,n;9 n = l) = 1 - V(p, n; 9 n = 0) = x n 6k ^ u ' ' ' ' Xn ~ l] . (4.19) 

The last relation follows from the trivial observation 

Cfe(^i) ■ ■ ■ , x n _i, x n ) = e^(xi, . . . , + x n e^_i(xi, . . . , x n ^i). (4.20) 

By the symmetry (of the symmetric functions), these probabilities extend to any 
orbital i : 

p(pn;9i = 0) = Z (P' n )U=° _ T l k=o e k(xi,---,x i ,...,x n ) ^ 

Z(p,n) J2k=0 e k( x l,---, X n) 

V(p,n;0 i = l) = 1 - P(p, n; t = 0)= x, E ^* (iCl ' / ' ' ' **' ' ' ' ' Xn) . (4.22) 

Herein, (xi, . . . , Xj, . . . , x n ) stands for the (n— l)-tuple obtained by removing x« from 
the n-tuple (x±, . . . , x n ). It is now clear that the average number of particles on the 
i-th orbital, denoted by 9 iy is just V(p,n; 9i = 1). In other words : 



X]fc=0 e k{ x li ■ ■ ■ 



9 t = ^k^*^y-> x *>---> x »> . (4.23) 

From (JHHJ) and ijOill this is also : 

e i = x i d s .QnZ(p,n)). (4.24) 
For p > n ()4.23|) gives the Fermi case 

H - T^. (4-25) 

which is consistent with (14.17(1 . So 9{ denotes the average number of particles on 
orbital i in the case of Fermi statistics. It is interesting to consider the deviation 
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when p = n — 1. We can express the average number of particles on orbital % in the 
case of A-superstatistics of order p = n — 1 by means of the Fermi averages d{ : 



ftf _ rp e f 

i-n; =1 ^ 

Clearly, these new averages are small deviations from the averages in the case of 
Fermi statistics. Also note that the deviation of the average on orbital i depends on 
the Fermi averages on all other orbitals. Similarly the average number of particles 
in the system for p = n — 1 is 

N(n, n)-n FT?, 
N( P = n- 1, n) = U n y ^ ( 4 - 27 ) 

where N(n,n) is the average particle number in the case of Fermi statistics. 
The average energy of the particles on the i-th orbital is given by 

Ei = eA = eiXids-OnZfan)), (4.28) 

and the average energy of the total system is 



n n 



E (P, n) = ^2 €i6i = e * x ^ ( ln Z (P' n )) 



1=1 1=1 

n p— 1 



f t j 

-— - y ejXj y £k{%ii ■ ■ ■ ■> ■ ■ ■ ■> %n)- (4.29) 
Z(p, n) ' ' — - 



i=l fc=0 



We can again express the average energy in the case p = n — 1 via the average energy 
in the Fermi case E(n,n), and the Fermi averages 9\ : 

E{p = n- 1, n) = 1 ' ^ fc=1 * llj=1 3 . (4.30) 

i-n7=i*/ 

Also other thermodynamical functions, such as the entropy S(p, n) and the heat 
capacity Cy(p,n), defined in terms of the grand partition function and the average 
energy, can be computed. The formulation of these expressions in terms of symmet- 
ric functions does not lead to further simplifications or insights, so we shall not deal 
with these. 



5 A-superstatistics in the degenerate case 

Let us consider as a particular example a Hamiltonian of the form 

n 

H = eJ2N t . (5.1) 

i=i 
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Thus we assume that all orbitals have the same energy, and let us furthermore 
assume that they also have the same chemical potential, i.e. \i\ — fi 2 — ■ ■ ■ — \i n — 
Therefore x\ — x 2 — ■ ■ ■ — x n — x, with 

x = exp(^^). (5.2) 
r 

In this case the orbitals label internal degrees of freedom of the particles such as 
spin, color, flavor, etc. The thermodynamical functions for this example follow from 
the formulas of the previous section, under the specialization x-i = x (i — 1, . . . , n). 
Since the number of terms in ek(x\, . . . , x n ) is given by (?), we have 



e k ( x,..^.,x) = y k J x ■ ( 5 - 3 ) 
n times 

Thus (gU) yields 

Z(p,n) = (5.4) 

For p < n this cannot be rewritten in a closed form; for p > n, this is simply (l + x) n , 
i.e. the GPF for a Fermi system with n distinct orbitals having the same energy. 
For p < n the sum in (15. 4|) can be rewritten as follows : 



k=0 V / k=p+l V 7 



x k 



1 + xT-l H i )x? +1 2 F 1 ( 1 > P ~ n n +1 ;-x), (5.5) 



kP + 1 / V P + 2 

where 2 F\ is the classical hypergeometric series 

^i( a '^)=£4#^ (d) k = d(d+l)...(d + k-l). (5.6) 

V C / k=Q \ C )k 

The first term in the right hand side of (j5.5J) is the Fermi GPF and therefore the 
second term gives the difference between Fermi statistics and A-superstatistics. In 
a sense, it describes the statistical interaction between the particles. 

Using Euler's transformation formula for hypergeometric functions, i.e. 



2 



F 1 l a ' b ;x)=(l-xr*- \fJ C a ' C b ;x), (5.7) 



C J \ c 

(15. 5|) can also be rewritten as 



Z(p, n) = (l + xr(l-( p n + J x^ 2 F X ( P + ^ + 1 ; -x) ) . (5.8) 
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So ()5.5|) gives the deviation from the Fermi GPF in additive form, and ()5.8|) gives 
the deviation in multiplicative form. 

The average number of particles follows from ()4.14|) : 



X 



k 



N(p,n)= ^ y fc . (5.9) 
Z^k=o \k) x 

Using the definition of hypergeometric functions, equation (|5.9jl can be rewritten as 

nx{ l + x) n-i _ (p + i) ( Mx^ 1 2 F X f 1 '^ 1 ; -x) 

N(p,n) = —J* i- (5.10) 

{l + xY-{; +l )x^ 2Fl (^--x) 

For p > n, this becomes 

— TlX 

iV(p >n,n) = - (5.11) 
The average number of particles on the i-th orbital follows from ()4.23|) : 

p-1 (n—X\ k 



(n-X\ 
,l^k=Q \ k I 



X 



°^ = x V A \ • (5-12) 
2^k=o \k) x 

For p > n, this sum simplifies to In the general case, simple properties of 

binomial coefficients lead to : 

- x fn — 1\ x p+1 

l+x \ p / (1 + x)Z{p, n) 

Obviously, this expression is the same for every orbital %. As a consequence, the 
average number of particles of the total system can be rewritten as 

nx fn — 1\ x p+l 



N(p } n) = -—-ni -. (5.14) 

l + x \ p J (1 + x)Z(p, n) 

It is interesting to consider an example. Let n = 5; we shall examine the depen- 
dence of the average number of particles on the i-th orbital 0j upon the variable 

y = — , where x = e~ v . (5.15) 

T 

In Figure 1, we plot 9i for p = 1,2,. ..,5. The case p = n = 5 yields the known 
Fermi-Dirac distribution function. For p > n, the distribution function is the same. 
For p < n = 5, the distribution function is different. The difference is most noticable 
for e < \i (or y < 0), as the average number of particles cannot exceed p/n. 
The case p = 1 and any n is also of interest : 

— 71 

N(l,n) = - — . (5.16) 
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N(l,n) is always smaller than 1, i.e. the system can accommodate at most one 
particle. When n = 1 this corresponds to the Fermi-Dirac distribution function. 
When n increases the average number of particles of the system increases also for 
fixed y = This description (p — 1, n > 1) corresponds to a system consisting of 
"hard-core bosons". Such particles appear (as mentioned in Sect. 3) in some models 
of condensed matter physics jTHj and nuclear physics [23J . This case is illustrated in 
Figure 2, where we take p = 1 fixed and let n vary. 



6 Equidistant energy levels 

Now we consider the Hamiltonian (J2.18j) with equidistant energies e«. Let the gap 
between the different energy levels be A > 0. Then 



ei = ei + (i - 1)A, (i = l,2, 



n) 



(6.1) 



We assume also that //i = //2 = • • • = = A*- Under these conditions the different 
orbitals correspond to different energy levels. According to notation of (|4.3|) . we 
have 

> - ej 



Xj = exp 



/i-ei\ / A 
exp | I exp I 



i-l 



xg 



i-l 



where we have used the notation 

' H - 6i 

x = Xi = exp ' 



(6.2) 



(6.3) 



Under this specialization the elementary symmetric functions simplify. For this 
purpose, consider their generating function (|4.7jl . Using [2U p. 26] one finds, 

(1 + xt)(l + qxt) • • • (1 + q n ^xt) = 

k(k-D/2 (<? n +1 ;<z)* 



^ g fc(fc-l 



)/2 



fc=0 



fc=0 



(?;?)/ 



(6.4) 



where [?] denotes the Gaussian polynomial [2U p. 26] and (a;q)k the g-raising 
factorials 25 

-ni _ (1 - g n )(l - g™" 1 ) • • • (1 - q n ~ k + l ) 
.hi ~ (l-g)(l-g 2 )---(l-g fc ) ' 
(a; q) k = (1 - a)(l - ga) • • • (1 - g fe_1 a). 
It follows from (jH?j) and (jHHJ) that 

-ni ,„ t , t . 1)/2 (g n ~ fc+1 ;g) 



e k (x, qx, q 2 x, • • • , g n x x) = g fe(fc 1)/2 



fc 



X A: = fc(ft-l)/2 



X . 



To write down the GPF, we use ()4.8|) . the specialization ()6.2|) and ()6.7|) 



Z(p,n) = ^g^~ 1 



)/2 



fc=0 



LA; 



x 



5> 

fc=0 



fc(fc-i)/ 



1)/2 (g- fc+1 ;g), fc 



(?; ?)/ 



X 



(6.5) 
(6.6) 

(6.7) 
(6.8) 



Using the formula 



-fc(fc-l)/2/ gU^k i-V n 'i 1)k 



one obtains 



Z[p,n) = ^ q ~J^\ -q n ) k x k 



k=0 



2^1 



q n ,q p 



-q x 



(6.9) 



(6.10) 



The last function is a terminating basic hypergeometric series |2*5] . 



a, 6 



£ 

fc=0 



{a;q)k(b;q) k x k 
(c;q) k (q;q) k 



For p = n ()6.8j) yields 

Z(n,n) = (-x; q) n 
The average number of particles in the system follows from ()4.14|) 



(6.11) 



(6.12) 



For p = n ()6.13|) becomes 



n-l 



N(n, n) = x — 



i=0 



(6.13) 



(6.14) 



The average number of particles on the ith orbital can be written in the form 

j-1 P- 1 k 



J(i-l)+(fc-0(fc-2-l)/2 



n 
k-l 



x 



(6.15) 



7 fc=0 /=o 

We used that : 

e k (x, qx, q 2 x, q l ~ 2 x, q l x, q^x) = ^(-1)V 

z=o 

(6.16) 

The conclusion from (j6.15|) is that the "population" of the orbitals depends mainly 
on their level i via q l ~ l with q = exp(— A/t) < 1 : as i increases, 9i decreases. 
Consider the case p = 1 and any n. Then 



n 
k-l 



x 



Z(l, n) = l + (l + q + q 2 + ... + q ^) e ^)h 
(l + q + ... + q n ~ 1 ) 



N(l,n) 



e (ei-M)/r _|_ (1 + q + _ _ _ + g n-l) 



(6.17) 
(6.18) 
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If q = exp(— A/r) <C 1, i.e., for large gaps between the energy levels or very low 
temperature, one can neglect all positive powers of q in (j6.18|) . What remains is the 
Fermi-Dirac distribution 

^' n ) W i^i - (6 ' 19) 
The expression for the average number of particles on orbital % reads 

q i-i 

9 i = —. — , i = l,...,n. (6.20) 

e (ei-M)/T + (1 + q + . . . + q n-i) ' y 1 

For very low temperatures, or a large energy gap A, ()6.20|) reduces to 

0i w and 6i w if i > 1. (6.21) 

Therefore if the system contains a particle, it is "sitting" permanently on the first 
(i.e. the lowest) energy orbital. This also explains why N(l,n) 9\. To illustrate 
these ideas, we plot the values of the average number of particles 9i on orbital i, in 
the case p — 1 and n = 5 (five equidistant energy levels), as a function of q, see 
Figure 3. For any value of q (0 < q < 1), 9\ > 9 2 > $3 > $4 > §5- For small values 
of q, 9\ is the large and the other averages close to zero. For increasing values of q, 
the averages on the other orbitals become larger. 



7 Concluding remarks 

In this paper we have studied the microscopic and thermal (macroscopic) properties 
of "free" particles, interacting only via statistical interaction. This interaction fol- 
lows from the Pauli principle of A-superstatistics : the system cannot accommodate 
more than p particles if the order of statistics is p, irrespective of the number of 
available orbitals (which may even be infinite). 

A-superstatistics is defined by the triple relations (j2.4|) - (j2.6j) which should hold 
for the creation and annihilation operators. The Fock spaces for A-superstatistics 
are naturally related to certain representations of the Lie superalgebra sl(l\n), and 
are labeled by a positive integer p referred to as the order of statistics. It is the 
mathematical structure of the Fock spaces that gives rise to the Pauli principle of 
A-superstatistics. 

It is shown that the creation and annihilation operators of A-superstatistics are 
fermion-like. Just like ordinary Fermi operators, the creation (resp. annihilation) op- 
erators anti-commute. However, they do not satisfy all traditional Fermi relations : 
only in the limit p — > 00 do the remaining relations tend to the ordinary Fermi 
relations. It may be unusual and unconventional to consider such alternative cre- 
ation and annihilation operators. However, we show that (representations of) such 
operators have already appeared in physical models, proving their applicability. 

In the second part of the paper we focuss on the macroscopic properties of A- 
superstatistics, for a free Hamiltonian. The usual thermal functions (grand partition 
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function, average number of particles, orbital distribution, average energy, . . .) are 
expressed in terms of elementary symmetric functions, and can be considered as 
deviations of the usual Fermi case when p < n. 

In addition to the general case, we have considered two specific examples. The 
case with identical energy levels per orbital (degenerate case) leads to further simpli- 
fications of the thermal functions, and the deviation from Fermi statistics (for p < n) 
becomes more apparent. The corresponding distribution functions are reminiscent 
of the Fermi-Dirac case, but there are also some striking differences (see Figures 1 
and 2). 

The case with equidistant energy levels is interesting from the mathematical 
point of view, since the thermal functions have simple expressions in terms of q- 
series (or basic hypergeometric series). This can be considered as a g-deformation of 
the degenerate case, and in fact the degenerate case can be deduced from the current 
one under the limit q — > 1. Also the situation with equidistant energy levels yields 
some interesting physical interpretations, e.g. concerning the orbital distribution 
(see Figure 3). 

In this paper, A-superstatistics was introduced mainly on the ground of mathe- 
matical considerations. Therefore, a natural next step would be to use this statistics 
for studying real physical objects or phenomena like an ideal A-gas or the quantum 
Hall effect, or statistics of cosmic rays etc. Such a study also includes a compari- 
son of the obtained results with existing experimental data or with the predictions 
of other noncanonical quantum statistics (parastatistics, fractional statistics, quon 
statistics, anyons statistics, Tsallis statistics or quantum group noncommutative 
statistics). We hope to return to this issue elsewhere. 
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Figure 1: Dependence of the average number of particles on the i-th orbital 9i upon 
the variable y = (e — /x)/r, for fixed n — 5, and p = 1,2,3,4,5. The distribution 
is independent of i since we consider the degenerate case here. For p = n = 5, 
the distribution function coincides with the Fermi-Dirac distribution; for p < n it is 
different and ^ cannot exceed p/n. 
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Figure 2: Dependence of the average number of particles N(p,n) upon the variable 
y — (e — /i)/r, for fixed p — 1, and n = 1, 2, 3, 4, 5, in the degenerate case. The graph 
of N(l, 1) is the closest to the y-axis, then 7V(1,2), etc. Observe that the graph of 
N(l, 1) coincides with the Fermi-Dirac distribution. 
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Figure 3: Graphs of the average number of particles 6i on the i-th orbital, for 
p — 1, n — 5, and % = 1,2,3,4,5. The distribution ^ is plotted as a function of 
q = exp(— A/r), where A is the gap between the equidistant energy levels. 
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